Multicritical Phenomena of Superconductivity and Antiferromagnetism in Organic 

Conductor k-(BEDT-TTF) 2 X 
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We study theoretically the multicritical phenomena of the superconductivity (SC) and antiferro- 
magnetism (AF) in organic conductors k-BEDT salts. The phase diagram and the experimental 
data on the NMR relaxation rate 1/Ti is analysed in terms of the renormalization group method. 
The bicritical phenomenon observed experimentally indicates the rotational symmetry, i.e., SO(5) 
symmetry, of the SC and the AF. The critical exponent x for the divergence of 1/Ti is well explained 
by x — u(z — 1 — 77) with the dynamical exponent z = 3/2 for the AF region while z — <p/v ~ 1.84 
at the bicritical point. These results strongly suggest that the origin of the SC is in common with 
that of the AF and that its symmetry is d-wave. 
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In strongly correlated electronic systems, it often hap- 
pens that superconductivity and magnetically ordered 
states touch or are near to each other. Heavy elec- 
tron superconductors WJ3|, and high-T c superconduc- 
tors H are examples of this phenomenon, where the 
Coulomb interactions are believed to be large com- 
pared with the bandwidth. In organic conductors, e.g., 
(TMTSF) 2 X (X:anion) @] and k-(BEDT-TTF) 2 X |||, 
on the other hand, the magnitude of the Coulomb in- 
teraction is smaller because the molecular orbitals are 
extended. However, the bandwidth W is also small 
(~ 1 eV) in these systems, and the electrons are half- 
filled when dimerization is considered. Therefore, it is 
expected that the Coulomb interactions play important 
roles in these compounds. In particular, in k-(BEDT- 
TTF)2X, recent experiments have revealed an interest- 
ing phase diagram including the bicritical phenomenon 
between the antiferromagnetism (AF) and the singlet su- 
perconductivity (SC) in the plane of temperature and a 
parameter p controlling the ratio of the Coulomb repul- 
sion U and the bandwidth W |Hq]- Above a character- 
istic temperature T*(U/W), the NMR relaxation rate, 
T^ , increases as the temperature is lowered in a man- 
ner independent of p. Below T*(U/W), on the other 
hand, T^ 1 diverges towards the AF transition tempera- 
ture, while it exhibits a spin-gap-like behavior on the SC 
side. 

It is a crucial theoretical issue to treat AF and SC 
in a unified fashion in the physics of strongly correlated 
systems. One of the most interesting proposals is based 
on the SO(5) symmetry H, i.e., the rotational symmetry 
in the 5-dimensional order parameter space of AF and 
SC. The critical phenomena are are useful for testing this 
idea because of their universality, and actually give the 
first evidence for it in these organic conductors. 

It is well known that the Ginzburg criterion T c| < 
(t) 2 specifies the critical region in 3D where the mean 
field theory breaks down. Here, a is the lattice constant 



and £ is the correlation length, whose ratio is determined 



roughly as | 



7^-) where A is the gap introduced 
in the electronic spectrum by the AF and/or the SC, 
and W e s is the reduced effective bandwidth by electron 
correlation. The mean field theory of the SDW for the 
AF and the BCS theory for the SC applies in the weak 
coupling regime, i.e., A <C W c ff, and the critical phe- 
nomenon is dominated by the mean-field-like behavior 
in this case. This occurs when the Coulomb interac- 
tion is weak (U <§C W) for the AF and/or the electron- 
phonon coupling is weak for the SC. In the strong cou- 
pling regime, on the other hand, we have a sufficiently 
large critical region and the critical behavior and the 
phase diagram are determined by the fluctuations beyond 
the mean-field theory. 

In this paper, we study theoretically the critical phe- 
nomena of the AF and the SC in k-(BEDT-TTF) 2 X from 
the viewpoint that this system is in the strong coupling 
regime and nontrivial critical phenomena are observed. 
The fluctuation of both the AF and the SC order param- 
eters are treated in terms of the renormalization group 
(RG) method, employing e = 4 — d expansion. The criti- 
cal phenomena occurs at finite temperature and hence is 
three dimensional. The anisotropy of the system is taken 
into account by proper rescaling and £ = y/£ x £y€z, an d 
the critical region is discussed later. Classification of the 
scaling trajectories leads to three types of multicritical 
phenomena (Fig. [l|): (a) tetracritical, (b) bicritical, and 
(c) tricritical phenomena. The bicritical phenomenon (b) 
is unstable towards tetra- and tricritical ones, and is re- 
alized only when there is SO (5) symmetry. Furthermore, 
the dynamic critical phenomena are studied, particularly 
for the NMR relaxation rate Tj -1 , and its critical expo- 
nent x is in reasonable agreement with the e-expansion 
result in both the AF region and at the bicritical point. 
Thus, the analysis of multicritical phenomena gives quan- 
titative evidence of SO (5) symmetry to a good accuracy 
without adjustable parameters, and it strongly suggests 



that the SC is realized by the Coulomb interaction as 
well as the AF and its pairing symmetry is d-wave. 

Let us consider a generic Ginzburg-Landau model for 
a system with competing AF and SC order: 



u,v,w up to order e are written as JO 



1 l -l 



H = d d r 

+ ulaf + 2w\a\ 2 \s\ 2 + v]^] , 



\\Vt\ 2 + \rM 2 + \\VA 2 



(1) 



where a and s are the order parameters of the SC and 
the AF, respectively. We normalize the order parame- 
ters so that the coefficients of the gradient terms are 1/2, 
and the other coefficients are roughly given in the unit 
where a = 1, W = 1 by r,, S %~^ c) , r± = ?"^ ) , 

u = ^j7 c , v = ^ F , while w represents the competition 
of the AF and the SC, and is expected to be positive. 
This model was proposed in refs. |lO|]ll| l in the context 
of the SO(5) theory for high-T c superconductors M, and 
was also used in ||12fl for the heavy fermion compound 
CezCu2Si2 (x ~ 1) to discuss the effect of disorder. 
This model has two types of mean field phase diagrams 
p3| . When the competition of the AF and the SC is not 
strong, i.e., w 2 < uv, there occurs the coexistence of the 
AF and SC, and the phase diagram looks like Fig.[y(a). 
When w 2 > uv, the AF and the SC are separated by a 
first-order phase transition line and the phase diagram 
looks like Fig. 0(b). One might be tempted to inter- 
pret the phase diagram and the bicritical phenomenon 
observed experimentally in terms of this mean field pic- 
ture. However, there are several experimental facts which 
imply that the system is in the strong coupling regime at 
least for the AF: (i) the large saturation value of the stag- 
gered magnetization M s = 0A/j,b Hi, which is nearly 
the value for the 2D Heisenberg AF, (n) fluctuation of the 
staggered moments are observed as the broadening of the 
NMR line shape well above T C! af an d (hi) the spin- 
gap behavior is observed in (TTi) -1 well above T Ci gc- 
Hence, the critical region near the bicritical point should 
be large provided that £af ~ 1- 
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FIG. 1. Schematic phase diagrams of the model (nj) as a 
function of t ~ (2rii + 3r±)/5 and g ~ (rii — r±). The thick 
and the thin lines represent first-order and second-order lines, 
respectively. "S.G." denotes a region with a spin-gap behav- 
ior. (A)-(C) correspond to materials discussed later in the 
text, (a) uv > w 2 : tetracritical, (b) uv = w 2 : bicritical, (c) 
uv < w 2 : tricritical. 



Let us proceed to the analysis of the fluctuation in 
terms of the RG method. The RG recursion relations for 
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It should be noted that the recursion relations for u, v, w 
do not involve ru , r±_ up to this order. There are six 
fixed points |15| , and a stable fixed point depends on the 
number of components n|| and n± of the vectors a and 
s. If n = nil + n± is smaller than n c = 4 — 2e + 0(e 2 ), 
a Heisenberg fixed point with u = v = w ^ is the only 
stable one. This fixed point becomes unstable when n 
exceeds n c , and the so-called biconical fixed point with 
unequal values of u,v,w becomes stable instead. 

In the present case with ri|| = 2, n± = 3, the only 
stable fixed point is the biconical one (ujg , i> B , u> B ) = 
27r 2 e(0.0905, 0.0847, 0.0536). However, not all points 
will flow to this fixed point through the RG (Fig||). 
There exists a curved surface, or a "separatrix" , 
F(u, v, w) — 0, which divides the entire parameter space 
into two regimes: one of convergence to the fixed point 
(u* ,v* ,w*), and the other in which the RG flows run 
away into an unstable region. Numerical analysis shows 
that this curved surface is well approximated by the sur- 
face uv = w 2 in the vicinity of the isotropic line u — v — 
w; thus, we can roughly say F(u, v, w) ~ uv — w 2 . 
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FIG. 2. Schematic diagram for the RG flow of the model 
(|l|) in the u/w-v/w plane. Broken lines represent the flow. 
Points "B" and "H" are the biconical and the Heisenberg 
fixed points, respectively. The shaded region is the unstable 
region of the model. 



If F(u,v,w) ~ uv 



> 0, the RG flow converges 



to the biconical fixed point (w B ,i> B , w B ) and the phase 
diagram shows a tetracritical behavior (Fig. Ha)), as is 
predicted by the mean-field approximation [ jl3[ . From 
the RG recursion relations of ru and r± with this bi- 
conical fixed point, we obtain j/ b = 0.5 + 0.132e and 
q b = — 0.0278e. The four second-order phase boundaries 
are written as \g\ ~ |t|^ B with the crossover exponent 
4>b = 1 + 0.135e, where t ~ r|1 5 r± and g ~ rn — r±_ are 
scaling fields corresponding to the temperature and the 
anisotropy between the AF and SC phases, respectively. 



We shift the definitions of t and g properly so that the 
tetracritical point corresponds to t = g = 0. The val- 
ues of the exponents are different from the ones in [jllj , 
which might be due to their approximation of smallncss 
of ; — -. A deviation of (u.v.w) from the biconical 

fixed point is irrelevant and affects only the size of the 
tetracritical region. 

If F(u,v,w) ~ uv ~ w 2 < 0, on the other hand, al- 
though the mean-field approximation predicts a bicritical 
behavior |13], it is not the case. Through the RG flow, 
u or v decreases rapidly to become negative, and the 
model (Of) becomes unstable. There are always higher- 
order terms which stabilize the system, though they are 
omitted in (pi). Accordingly, the phase transition between 
the disordered and the ordered phases becomes first order 
(fluctuation-induced first-order transition [pi), but if the 
quadratic anisotropy \g\ is sufficiently large, the transi- 
tion becomes second order again |l7J]. Thus, a first-order 
transition line between two ordered phases branches at 
the triple point and extends until these two branches ter- 
minate at tricritical points (Fig. 0(c)). If (u, v, w) is near 
the surface uv — w 2 , the length of the first-order lines 
on the normal-SC or the normal-AF phase boundaries 
will be small, since it takes a long time to flow into an 
unstable region. 

Only in the case uv — w 2 does the RG analysis pre- 



dict a bicritical behavior (Fig. Ill(b)) . This behavior is 
governed by the Heisenberg fixed point ujj = t;Jj = w^ = 



2n e/13, which is stable only on the surface uv 



and 



the corresponding exponents are v^ — 0.5 + 0.135e, 0h = 
I + 0.192e,a H = -0.0385e. The experimentally ob- 
served bicritical behavior, combined with the above dis- 
cussion, strongly suggests uv = w 2 to a good accuracy, 
which corresponds to the rotational symmetry in the five- 
dimensional order parameter space of (a, s). It is inter- 
esting to note that only if uv = w 2 is the model (|I|) 
smoothly related to the SO(5) NLcr model in pjlq|, in 
the limit —r» , — r±, u, v, w — > oo with their ratios fixed. 

Now let us consider the dynamic critical phenom- 
ena Q. The NMR relaxation rate ||§| exhibits a di- 
vergence toward the transition in the AF side, or a spin- 
gap behavior in the SC side, while these two cases show 
a similar behavior above T*. This is reminiscent of the 
bicritical crossover behavior of the model (pi) . The NMR 
linewidth is proportional to (T — T c ^-p)~ x on the AF 
side of the normal phase with x — v{z — I — 77), where 
z is a dynamic critical exponent [ ^9[ . When the system 
is in the vicinity of T c .af, but not very near the bicrit- 
ical point, its dynamic critical behavior is governed by 
that of the AF isotropic Heisenberg model. Macroscopic 
varibles describing slow dynamics are the staggered mag- 
netization s and the conserved uniform magnetization to, 
giving z — d/2 = 3/2. Thus, the exponent is xaf = 0.315 
up to 0(e 3 ) (g0[. On the other hand, when the bicritical 
fixed point governs the dynamic critical phenomena, the 
value of z changes, and so does the exponent x. In the 
bicritical region, the SC order parameters a enter the set 



of macroscopic variables. The free energy F has a scaling 
form: 



F(t,g,a, s,m) 



t 2 - a <f> 



g a s m 



(5) 



where — 2 — a — 0h- t is proportional to T — T Ci bp, 
where T Ci bp denotes the temperature at the bicritical 
point. Following the power-counting argument of the bi- 
critical behavior in the spin- flop AF pij| , we get z = 
^h/^h ~ 1.84, resulting in iebp = 0.584, where we used 
4> H - 1.313, v H - 0.714 up to 0(e 2 ) @. This value of 
xbp at the bicritical point increases as a function of n: 
x BP = 0.528, 0.558, 0.584, 0.607 for n = 3,4,5,6, respec- 
tively. On the other hand, in approaching the SC phase 
it does not diverge but exhibits a spin-gap behavior be- 
low a characteristic temperature ig G ~ Igl 1 ^ 11 because 
of the singlet formation. Throughout the entire criti- 
cal region of the normal phase, T^ 1 has a scaling form 



t; 



t~ XBP f(g/t 4 ' H ). 
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FIG. 3. Log-log plot of T-f 1 vs (T - T c )/T for 

(A) k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 (solid squares), and (B) 
deuterated K>(BEDT-TTF) 2 Cu[N(CN) 2 ]Br (open squares). 
(Data from M.) 



Fig. shows the log-log plot of T l l vs 



T-T c 



for 

(A) k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 (solid squares), and 

(B) deuterated K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br (open 
squares). The former is located in the AF region slightly 
distant from the bicritical point while the latter is nearly 
at the bicritical point, as shown in Fig.m(b). We selected 



T-T a 



instead of 



T-T c 



as the abscissa, because the crit- 



ical temperatures of these compounds are much lower 
than their classical values (~ J ~ 500 K), indicating 
that they are close to a quantum critical point. This 
scaling by ~Z , near the quantum critical point is ob- 
tained from an RG study of a quantum SO (5) NLcr model 
with a quadratic anisotropy, and a detailed discussion 
will be given in a subsequent paper [E2J. The proximity 
to the quantum bicritical point is also consistent with 
the strong correlation, because T c is suppressed by the 
quantum fluctuation driven by the Coulomb interaction. 



The estimated slope, i.e., the critical exponent x, in the 
region 1 K < T - T CjA f ^ 10 K gives x AF ~ 0.30 ± 0.04 
for (A) and xbp ~ 0.56±0.04 for (B). The values of cc are 
in reasonably good agreement with the above theoretical 
values (xaf = 0.314, xbp = 0.584), which supports that 
the critical region is wide (~ 10 K) and £af ~ 1< Several 
points should be noted, (i) The fitting with T T r Ta under- 
estimates the values of x as x ~ 0.24 for (A) and x ~ 0.40 
for (B); thus, it guarantees the validity of our fitting with 
^=±k (ii) Note that the data in T - T c < 1 K seems to 
deviate from the above fitting. We discard these data 
in the fitting for the following reasons. First, T c cannot 
be determined more accurately than AT C ~ 0.5 K due to 
experimental limitations. This uncertainty largely affects 
the fitting within T — T c < 1 K. Second, in real systems, 
SO (5) symmetry is approximate, and there is always a 
small deviation from uv — w 2 = 0, which alters the crit- 
ical behavior very close to T c . Here, it is important to 
note that the increase of the deviation through the RG 
flow is slow; in the linearization of (g)-(Q), the exponent 
e/13 of the growth is much smaller than the magnitude 
of the other two exponents — 8e/13, — e. Therefore, the 
bicritical behavior governed by the SO (5) fixed point can 



be observed for a relatively wide range of \uv 



w 



For 



example, 



< 0.1 • max(wu,u> 2 ) is sufficient to 



explain the observed bicritical phenomenon. 

Let us consider whether it is possible to have such a 
large (~ 10 K) critical region in the quasi-2D system. 
Following [E3| , we can evaluate the critical region of the 
quasi-2D Heisenberg system as 



\T-T C 
T, 



< 



kT r 



1 



j\\ hj\\/j±v 



(6) 



where the in-plane exchange Ju is larger than the inter- 
plane exchange Jj_. Thus there is only logarithmic de- 
pendence on J±/Ju, and it does not significantly reduce 
the width of the critical region. In k-(BEDT-TTF) 2 X in 
a metallic region, the anisotropy of the conductivity is 
fry/o-j. - 100, implying t\\/t± ~ VlOO and J\\/J± ~ 100. 
Thus, the width (g) is not extremely small (~ 0.2 — 0.3), 
and in the AF region (T c ~ 30K) the observed size of the 
critical region 5 K — 10 K is reasonable. In contrast, in 
quasi-lD systems, the width in (m is J±/J\\; the quasi- 
one-dimensionality is effective in reducing the width of 
the critical region. 

On the SC side, the experimental estimation of the 
coherence length £sc is still controversial, e.g., for the 
SC compound k-(BEDT-TTF) 2 Cu(NCS) 2 , the calcu- 
lated in-plane coherence length £sc ranges from 31 A 
p4f to ~ 180 A J25|. However, the rotational symmetry 
between the AF and SC further suggests that the cor- 
relation length £ sc is also short (~ 1). This is because 
when £sc ^ £af, the scaling trajectory in Fig. || starts 
from the initial point (u <C v) far away from the Heisen- 
berg fixed point, and the deviation from the separatrix 
is easily magnified as the length scale becomes larger. It 
is also inferred that the observed spin-gap behavior is 



due to the large fluctuation of the SC order parameter 
in the SC compound k-(BEDT-TTF) 2 Cu(NCS) 2 ((C) in 
Fig. |(b)). 

One might wonder that the easy-axis anisotropy re- 
duces the number of the components of the AF order pa- 
rameters and stabilizes the bicritical phenomenon. How- 
ever, the spin-orbit interaction is negligible in organic 
systems, and the spin anisotropy energy due to dipole- 
dipole interactions is estimated by the AF resonance |29] 
to be - 10~ 4 K for (TMTSF) 2 X. Assuming a similar 
value for k-(BEDT-TTF) 2 X, and that £ A f ~ 1, this 
anisotropy is negligible in the temperature region of in- 
terest. 

Finally, the implications of these results are discussed 
below. The rotational symmetry between AF and SC 
order parameters suggest that the mechanism of the 
AF and SC is common J27) and that the underly- 
ing microscopic quantum model has enhanced dynam- 
ical symmetry, i.e., SO (5) ||. Therefore, it is likely 
that the symmetry is d-wave. We believe that the (a) 
half-filling, (b) intermediate Coulomb interaction, and 
(c) nearly two-dimensional Fermi surface in k-(BEDT- 
TTF) 2 X makes the SO(5) symmetric model a promising 
candidate. When these conditions are largely violated, 
e.g., by doping carriers and/or the application of an ex- 
ternal magnetic field, the SO (5) symmetry is broken and 
the bicritical phenomenon turns into tetra- or tricritical 
behavior 18 1. 
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